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The magnetic field responses of the zero-bias conductance peak (ZBCP) in tunneling spectra 
of high Tc cuprate junctions are studied theoretically. Our calculation is based on the lattice 
Green's function method and takes the realistic electronic structure of the high Tc cuprate 
into account. In marked contrast to previous works, it is shown that the critical magnetic 
field strength He exists for the splittings of the ZBCP's to be discernible. He is almost 
proportional to the product of the magnitude of pair potential, transmissivity of the junction, 
and the inverse of the Fermi velocity parallel to the interface ( l/vpy )■ The calculated -ffc's 
for the hole-doped superconductors are higher than those for electron-doped ones because of 
relatively large magnitude of pair potential and the small magnitude of vpy originating from 
peculiar shape of the Fermi surface. 
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1. Introduction 

Nowadays, almost all of the high-Tc cuprates are clar- 
ified to be c?-wave superconductors^"^-' and several exper- 
iments reported the existence of phase coherent phenom- 
ena peculiar to the d-wave symmetry. One of this is so 
called the zero bias conductance peak (ZBCP) observed 
in the tunneling spectroscopy in high Tc cuprate junc- 
tions. It was clarified that ZBCP is a direct consequence 
of dx2_y2 symmetry of the pair potential. ^"'^'^^ The pres- 
ence of ZBCP is a manifestation of the formation of the 
Andreev bound state (ABS) at the Fermi energy (zero- 
energy) near a specularly reflecting surface^^ or inter- 
face, when the angle between the lobe direction of the 
dx2_y-2--wsLYe pair potential and the normal to the inter- 
face is nonzero. This state originates from the interfer- 
ence effect between the injected and reflected quasiparti- 
cles at the surface or interface and the sign change of the 
d^2_j^2-wave pair potential. There is a large number of 
studies on ZBCP,''' ^^"^^^ and its related problems. ^^^^^^ 
However, the experimental results of the dependence of 
ZBCP on an external magnetic field are not converged 
yet. The splitting of the ZBCP in a magnetic field has 
been reported in some experiments"'^^' ^^"^''^ while oth- 
ers can not reproduce it.^^"^*-' We can consider whether 
splitting of ZBCP is observed or not is quite sensitive to 
the experimental situations and doping concentration of 
actual high Tc cuprates. 

The effect of screening currents which is so called 
Doppler shift of the energy levels of quasiparticles is con- 
sidered to be the most promising idea to understand the 
origin of the peak splitting.^^'^^'^^^ It was pointed out 
that in an applied magnetic field H, screening currents 
shift the ABS spectrum^^^^^^ and lead to a splitting of 
ZBCP that is hnear in H at low fields.^^) The previous 



theory^^^ predicts the splitting of ZBCP in any case and 
is consistent with some experimental results. ^^"^^-' How- 
ever, the absence of ZBCP splitting has not been clari- 
fied yet. We must need further study which takes into 
account the electronic structures peculiar to actual junc- 
tions of high Tc cuprates. Recently, four of the aurhors 
study about the effect of transmissivity of the junction on 
the sphtting of ZBCP.32) It was found there is a critical 
value of the magnetic field He above which ZBCP splits 
by external magnetic field below which it does not.^^^ 
The critical field is almost proportional to the width of 
ZBCP in the absence of magnetic field and equivalently 
to the transparency of the junctions. However, in this 
paper, the shape of Fermi surface is assumed to be cylin- 
drical and the reduction of the magnitude of the Fermi 
velocity parallel to the interface in the actual Fermi sur- 
face is not taken into account. Since the degree of the 
Doppler shift is determined by the Fermi velocity parallel 
to the interface, it is natural to expect that the shape of 
Fermi surface influence significantly on the ZBCP split- 
ting. In order to understand the actual line shape of the 
tunneling conductance, it is indispensable to calculate 
tunneling conductance taking into account the effect of 
the shape of Fermi surface. 

The aim of this paper is to predict the tunneling con- 
ductance of normal metal - high Tc superconductor junc- 
tions in the presence of magnetic field taking into account 
the shape of Fermi surface. In the present paper, based 
on the Kubo formula, '^^"'^^•' we will clarify that there is 
a critical value of the magnetic field -ffp, which crucially 
determines the peak structure of the tunneling conduc- 
tance. For H > He, ZBCP splits into two by the mag- 
netic field.26) On the other hand, for H < He, ZBCP 
does not split at least explicitly and only the height of 
ZBCP is reduced. It is revealed that He is proportional 
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to the product of the width T of ZBCP ioi H ^ and 
the inverse of the magnitude of Fermi velocity in c?-wave 
superconductor parallel to the interface (vpy)- Since the 
magnitude of T is roughly proportional to the product of 
the transmissivity of the junctions and the magnitude of 
the pair potential, He is enhanced (suppressed) for high 
(low) transmissive junctions and large (small) magnitude 
of the pair potential of c?-wave superconductors. We also 
clarify the asymmetric behavior of He with respect to 
the carrier type of the high Tc cuprates mainly due to 
the difference in the shapes of the Fermi surface. 



the applied magnetic field. The quantity ky is replaced 
with ky — eHX/H. In the following, the applied magnetic 
field is normalized by Hq = 0o/(7''^Co'^) with 0o = h/{2e) 
and ^0 = fT'Vpo/{AoQTr). The order of the magnitude of 
Hq becomes about 1 Tesla when we choose and A 
as lOAand 1500A, respectively, for hole-doped cuprates 
with S = 0.1, where Aqo = 0.12to and vfq = O-lGtoa' /h 
with the lattice constant a' of cuprates and the unit to 
of energy. 

As for the normal region (N and /) , we use the single- 
orbital tight-binding model. The Hamiltonian is 
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Fig. 1. Schematic illustration of the geometry of the normal metal 
(N)/ insulator (I)/ d^2_y2-wave superconductor (D) junction. N 
and I are stacked along (100) axis of a square lattice while D is 
stacked along (110) axis. Lattice spacing is a for N and I and 
a'(= a/V2) for D. 



(3) 

l,m,(T 

where t is the hopping integral between nearest neigh- 
bor sites and vi^m is the on-site potential at a site (/, m). 
Since the system considered has the translational invari- 
ance in y-direction, the conductance Ts{eV) of the junc- 
tion is given by Kubo formula'^^"'^'^-' using the Green's 
function G{E) as 
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2. Formulation 

We consider a N/ 1 /D junction on 2-dimensional lat- 
tice system, where N and / are stacked along (100) axis 
of square lattice with lattice spacing a and I is connected 
to the (110) surface of D (see Fig. 1). In order to describe 
the d-wave superconductor, we use the t-J model which is 
one of the promising models for high Tc cuprates. By us- 
ing a mean-field approximation based on the Gutzwiller 
approximation, the model can be mapped into a BCS 
like mean- field Hamiltonian^s) + A with 
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where summations over (• ■ ■ ) and (••■)' run nearest- 
neighbor and next-nearest-neighbor pairs, respectively. 
The hopping parameters ii and ^2 include so-called 
Gutzwiller factor and the pair potential Ajj has d^^i-y^ 
symmetry. In the middle of the superconductor, Ajj 
is Ao(<5) for j = i ± (1/2,1/2), -Ao((5) for j = i ± 
(1/2,-1/2), and zero for others, where 5 denotes the 
doping concentration of hole or electron. Due to the 
translational invariance along the y direction, the cre- 



ation and annihilation operator c. 



and Cj ^ can be 



expressed as cj ^ = Cj^^a-iky) ex.p{—ikyjya). Since 
the coherence length of the high Tc cuprates ^ is much 
smaller than the penetration depth of the magnetic field 
A, and the effect of the magnetic field on the normal 
metal is not important, we can choose the spatial depen- 
dence of vector potential as Ay{x) = —HX, where H is 



ky ky ^ ky ^ ky 

^1,1 ^1+1,1 + 1 ^ '-'1 + 1,1+1'-' 1,1 



G = G{E - iO) - G{E + iO) , 



(4) 
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where the relation ci,a-{ky) — J^m'^i,' 



(6) 

a- exp{ikyma) was 



used and G^ is a 2 x 2 matrix. First, we calculate the 
isolated Green's function of the d-wave superconductor 
and normal metal. Then we obtain the combined Green's 
function by solving Dyson's equation iteratively and fi- 
nally we obtain the tunneling conductance Ts{eV).^^'^'^'' 
We also calculate the corresponding quantity TNieV), 
where the d-wave superconductor is in the normal state 
and the pair potential is, then, zero. Hereafter, we will 
look at the normalized value arieV) = rs{eV)/Tj\[{eV). 

In the actual calculation, we choose material param- 
eters in the d-wave superconductor determined by ap- 
plying the Gutzwiller approximations to t-J model. For 
both the normal metal and insulator, we use t — lOto 
and 11 = 0. As regards w;,m, we choose vi^m = 20to for 
normal metal, while for insulator we choose w/.„i = 35io 
for the junction with high transmissivity, and vi,m = 65to 
for junctions with low transmissivity. The corresponding 
transmissivities 7 are given by 0.4 and 0.1, respectively. 
The thickness of the insulator is chosen as one atomic 
layer. 

3. Numerical calculation and results 

First, we study the electron over-doped case, i.e., 
5 = 0.2, (Fig. 2), where = 0.4150to, ^2 = 0.1333io, 
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Fig. 2. Tunnel conductance crx{eV) of N/ I / D(clectron dope, 
S = 0.2) junctions with (a) high (7 = 0.4) and (b) low (7 = 0.1) 
transmissivities. Results obtained for H = 0.0, 0.3, 0.5, and 1.0 
in a unit of Ho axe plotted by dashed, dotted, chained, and solid 
lines, respectively. 



Fig. 3. Tunnel conductance axieV) of N/ I / D(hole dope, S = 
0.2) junctions with (a) high (7 = 0.4) and (b) low (7 = 0.1) 
transmissivities. Results obtained for H = 0.0, 0.3, 0.5, and 1.0 
in a unit of Ho are plotted by dashed, dotted, chained, and solid 
lines, respectively. 



/i = 0.0214^0 for the d-wavc superconductor.'^*^ The 
ZBCP sphts into two as the magnetic field increases (see 
Fig. 2). The degree of the splitting is crucially influenced 
by the transmissivity of the junctions. With an increase 
in transmissivity of the junctions, the critical value of 
magnetic field He above which ZBCP splits into two 
decreases. 

In Fig. 3, corresponding plot for hole over-doped case 
{S = 0.2) is shown, where material parameters are chosen 
as <i = 0.4117to, t2 = -0.1333*0, and /x = -0.4601io-^^^ 
For hole-doped materials, the degree of ZBCP splitting is 
weakened, and ZBCP does not split into two by magnetic 
field except for low transmissivity and high magnetic field 
junctions. As compared to Fig. 2, the magnitude of He 
is drastically enhanced. 

To clarify this point, 8 dependence of He = Hc{5) is 
plotted in Fig. 4 both for electron and hole-doped cases 
with low and high transmissivity of the junctions. The 
magnitude of Hc{5) is a decreasing function with 5 both 
for electron and hole-doped cases. In the following, we 
will consider a microscopic origin of Hc{S). 

In the presence of magnetic field, the energy of the 
quasiparticle E is substituted with 



given by 



E — evpyAyix) = E + 



HqVfQ 



(7) 



Since positive H is chosen, quasiparticles energy E in- 
creases (decreases) by H for vpy > (vpy < 0). The am- 
plitude of this energy shift is HAq < \ vpy \>/{Hovfo), 
where <| vpy \> denotes the average of the y compo- 
nent of the Fermi velocity on the Fermi surface. By 
comparing the energy shift and the width F of ZBCP 
with <tt(^) = ^ariO) for H = 0, we can analytically 
estimate critical value He which we call Hth{S) to dis- 
tinguish from Hc{5) in Fig. 4. The resulting Hth{S) is 



HthiS) 



CTHqvfo 

2Aoo <\vFy\> 



(8) 



where C is a constant. To confirm the correspondence 
between HdS) and Hth{S), we plot Hc{6)/Hc (0.10; 
electron-dope) and Hth{S)/Hth (0.10; electron-dope) in 
Fig. 5 for low and high transmissivity of the junctions 
with both electron and hole-doped cases. 

The coincidence between these two values are fairly 
well, and we can expect Hc{6) ~ HthiS). As discussed 
above, there are two crucial factors which determine the 
magnitude of He- First one is the magnitude of F which 
is almost proportional to the product of the transmissiv- 
ity of the junctions and the magnitude of the pair poten- 
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Fig. 4. i?c'((5)/-ffo ca-lculatcd for junctions with high (7 = 0.4) and 
low (7 = 0.1) transmissivities as a function of the doping rates 5 
of hole (triangles) and electron (circles). Solid and dashed curves 
are guides for eyes. Ho is about IT. 
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Fig. 5. Normalized He and Hth calculated for junctions 
with high (7 = 0.4) and low {7 = 0.1) trans- 
missivities as a function of the doping rates <5 of (a) 
hole and (b) electron. ifc? (i5)/if £7(0.1; electron-dope) and 
^th(<5)/^^th(0.1; electron-dope) are plotted by solid and open 
symbols, respectively. Solid and dashed curves are guides for 
eyes. 



tial Ao((5). The magnitude of Ao((5) is enhanced with the 
decrease of 5 as seen from Fig. 6. The second one is the 
magnitude of <| Vpy \>, which is reduced significantly 
for hole-doped high Tc cnpratcs (sec Fig. 6). In the light 
of our calculation, the absence of ZBCP splitting in the 
tunneling experiment of hole-doped high Tc cuprates is 
fully plausible, since the magnitude of He is enhanced. 




in detail where the realistic Fermi surfaces of high Tc 
cuprates are taken into account. We concentrated on 
the splitting of the ZBCP due to the magnetic field for 
(llO)-oriented junctions. We can define a critical mag- 
netic field He, where only for H > He, ZBCP splitting 
occurs by magnetic field. It is revealed that He is 
proportional to the product of the transmissivity of the 
junction,^^) the magnitude of pair potential, and the in- 
verse of the averaged Fermi velocity parallel to the in- 
terface. In the light of our theory, the absence of the 
splitting of ZBCP reported in several experiments^''"^^) 
is due to the larger magnitude of He originating from 
the low doping concentration or high transmissivity of 
the junctions. Systematic studies of tunneling experi- 
ments changing doping concentration in the presence of 
magnetic field are desired. 

In the present paper, the effect of the impurity scat- 
tering or surface roughness is not taken into account. 
These effects broaden the ZBCP=^9'^2,43) ^^^^ ^^le magni- 
tude of Hc{5) is expected to be enhanced. It is actually 
an interesting problem to clarify the influence of the im- 
purity scattering on He {5). 

In the present paper, spatial dependence of the pair 
potential is not taken into account. It is revealed that if 
the symmetry of the pair potential is pure d-wavc;. the 
low bias property of the tunneling conductance is not 
influenced by the spatial modulation of the pair poten- 
tial.*'' However, in the presence of the broken time re- 
versal symmetry state (BTRSS) near the interface, where 
the mixing of subdominant s- or dxy-wave component as 
the imaginary part of the pair potential to the predom- 
inant da;2_y2-wave component^^''^^"*^'*'*'*^) occurs, the 
ZBCP splits into two even without magnetic fleld. It is 
also one of interesting problems to clarify the magnetic 
field dependence of the conductance in the presence of 
the BTRSS.46) 

Finally, we must comment about recent experiments 
by Dagan and Deutsher.^'^'^'') They have found that the 
width of the splitting of ZBCP is thickness independent 
for films thinner than the penetration depth. This im- 
plies that the magnetic field splitting is not primarily 
a Doppler shift effect due to the Meissner screening cur- 
rents. We must seek for other possibilities including field 
induced modification of the pair potential.^^'^*'''''^ 
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lutional Science and Technology (CREST) of the Japan 
Science and Technology Corporation (JST). J.I. acknowl- 
edges support by the NEDO project NAME. The com- 
putational aspect of this work has been performed at 
the facilities of the Supercomputer Center, Institute for 
Solid State Physics, University of Tokyo and the Com- 
puter Center. 



Fig. 6. The magnitude of vpy and Ao(5) as a function of 5. 

4. Conclusion 

In this paper, influence of the magnetic field on the 
tunneling conductance in N/I/D junctions is studied 
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